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A non-linear Penner type interaction is introduced and studied in the random matrix model of 
homo-RNA. The asymptotics in length of the partition function is discussed for small and large 
j ' (size of matrix). The interaction doubles the coupling (w) between the bases and the dependence of 

I ' the combinatoric factor on (v, N) is found. For small N, the effect of interaction changes the power 

^ , law exponents for the secondary and tertiary structures. The specific heat shows different analytical 

' behavior in the two regions of A'', with a peculiar double peak in its second derivative for A^ = 1 at 

low temperature. Tapping the model indicates the presence of multiple solutions. 

S ■ 

^ . PACS numbers: 02.10.Yn, 02.70.Rr, 05.40.-a, 87.10.-e 

I In the fundamental understanding of RNA folding combinatorics, the study of exact enumeration of RNA secondary 

^^^^ ' structures with crossings (pseudoknots) is an important ongoing research direction [ij. In this context, Orland 

I— 1 1 and Zee proposed a random matrix- field theoretic modelM which addressed the problem of exact RNA structure 

i-C ■ combinatorics which with certain simplifying assumptions [1[ enumerated all possibleplanar and non-planar structures. 

^ _ In this model, the idea of introducing external linear interaction was explored in with the objective to observe 

' structural and other statistical changes (if any) . Even with this simple interaction, it was shown that the interaction 

I ' imposed additional constraints on the RNA chain which could physically imply changes in temperature, applied 
pressure, proximity with ions @. So, within this matrix model framework, an important interaction to study is the 

> [ Penner type that may capture interesting properties such as: effects of interactions with more complex molecules 

^ • • and biomolecules, multiple solutions, frozen (glassy) states The Penner matrix models appear in the context of 

d [ disordered systems Q (initially to calculate correct fluctuations for the conductance), string theory Q (to get accurate 

^ • critical exponents for quantum gravity) and spin glasses (as mappings to high temperature p-spin glasses) . In this 

I ■ letter studying the Penner type interactions increases the moduli space of structures from Mg^ to Mg.n where g is 

' ^ 1 the genus of the surface and n gives the number of faces or punctures of the Riemann surfaces. Thus the generalized 

g partition function of RNA matrix model (of length L: — 1, L) with a logarithmic interaction in the action is 

Q - given by 



O 

in 
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where (/)i's are L independent {N x N) random (symmetric) hermitian matrices placed at each base site in the chain 
with the interactions contained in Vij (see for all other notations and conventions) . The normalization constant 
is given by 

O ■ 

>< 

' and jj'TrY\^{l + ipi) is the characteristic observable of the model. With the simpliflcations, Vij = v and Wi 



and a series of Hubbard Stratonovich Transformations |4l, Eq. (1) becomes 



(3) 



where cr is an (N x N) matrix, the potential is given by V{(t) — [§cr^ — ^loga'^] (with /i = t" — 1/w) and 
the normalization is 

MN) = j dae-^^^^^^"-^'°s^"l (4) 
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These potentials are of the Gaussian Penner matrix models ^ and will be solved along those lines. The spectral 
density pn{z) of the matrix a at finite N is 



Pn{z) 



N 



(5) 



where z are the eigenvalues of a. Defining G{t,N) — X^lLo ^l{N)jj as the exponential generating function of the 
partition function [s!, and using the identity dzpiq{z) — 1 gives 



+ 00 



G{t,N)^ dzpN{z)exp"^^+^\ (6) 

J — OO 

To solve G{t^N), the expression for spectral density is found using the orthogonal polynomial method (Deo, [9]). 
For these models, the orthogonal polynomials are given by Vn{z) = z^ + /.o. which satisfy the orthogonality condition. 
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dze-^'"^'^Vn{z)V,n{z) = K5„ 
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For the symmetric Gaussian Penner matrix model, orthogonal polynomials split into even and odd polynomials. 
The even set obeys the orthogonality condition 



where y = N' ^ N/2, vn{y) = 2Vo(z) = py + ■■■■ and t' = (t" - j^). The odd ones obey 



(8) 



dye-'''^''°y-'''°^y¥n(.y)P,n{y) - hnSnrn, 
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where i' = {t" + 5^). It is sufficient to work with either of the two polynomials as each one can completely 
determine the recursion coefficients independent of the other. The normalized even orthogonal polynomials from 
solving the orthogonality relations are 



n\{N'p) 
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which in the case of odd polynomials has t' instead of t'. The kernel for this function is defined by K{yi,yj) — 
'l2n=o '4'2n{yi)ip2n{yj) which gives the normalized spectral density for the even polynomials as p%{y) — j^K{y, y) (the 
superscript e represents even). Thus the spectral density for large N limit is 
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where y — z"^ \s considered. Using the relation p%{z) — zp'^j^{z^) (Tan; Deo in Q) and substituting p^(z^) in the 
generating function for the even polynomials, Eq. (6) gives 



G^it,N) 
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To solve this integral, a simplification is made. N't' — 1/2 (—1/2 for the odd) which requires t" ~ 2/N {—2/N 
for the odd). Substituting this and using the relation _ff2n+i(a;) = {—l)"'2'^"-^^nlxLl/^{x^) (for odd polynomials, 
H2n{x) = (-l)"22"n!L,;^/2(a;2)) in Eq. (12) gives 
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TABLE I: The Table lists even partition functions for L upto 6 for the matrix model with logarithmic interaction. 



L Zl{N) 

1 1 

2 l + ^(2+i) 

3 l + 3w(2+i) 

4 l + 6«(2+i) + 2i;2(4+f) + ^ 

5 l + 10t;(2 + ^) + 10^;2(4+|p) + f^ 

6 1 + 15«(2 + i) + ■iQv\A + A) + 5„3(8 + ^ + ^) 
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Using the formula [Hk{x)f - ^f^^ (^!)!|l±iy2i(x) and the integral da;e-(--«)'i/„(x) = Vny"2", G^(t,iV) 



becomes 



= 2 ,14, 

n=0 /=0 ^ ^ ^ 

which differs from odd polynomials in X]f"o (^i")- '^^^ partition functions can be obtained from Eq. (13) (Table 
I). The general form of the partition function for any N't' requires a rigorous mathematical analysis and is left for 
a future work. In this model, A'' plays a dual role of contributing to the strength of the external interaction and a 
genus identification parameter of the structures. However with N't' = (1/2, —1/2), the genus characterization cannot 
be extracted systematically. In the model of 0, Q , the structures (and their genus characterization) were obtained 
on the moduli space of a zero puncture Riemann surface A^g,o- Constructing the matrix model of RNA with a 
logarithmic interaction generalizes the study of RNA structures to that of n-punctured Riemann surfaces, M.g,n- The 
Euler characteristics for these models are given by (T^ — i? + n) = (2 — 2g) where V and E give the number of vertices 
and edges respectively and includes the additional factor of faces or punctures n. The genus characterization of the 
structures obtained from the model is therefore changed from Q . For this model the Feynman diagrams are given by 
the fat-graphs [Si]. 



ASYMPTOTIC ANALYSIS OF THE PARTITION FUNCTION 

The asymptotic behavior of the partition function at large length is found numerically as in jBj. The analysis is 
divided into two parts, (I) the estimation of structure combinatoric factor and (II) the determination of (secondary 
and tertiary) power law exponents with L. 



Combinatorics 



The combinatorics of the structures is given hy X = {2^J f{v, N) + 1)'" where (i) X = 3^ for the model in [3] and (ii) 
X = {3 — a)^ in (here f{v, N) = w = 1 in both the cases). The analysis is done for lengths upto L = (40, 80, 160). 
For each length, different v (1 to 6) are considered and for each v different values of N (1 to 100000) are chosen so 
as to observe the effect of interaction. In order to determine the form of f{v,N) for this model, /n[Z£(A^)] is plotted 
with L for different {L,v,N). The linearly fitted slopes and intercepts are found. For a given v, 

(i) /(v, A^)'s for each N are found from the slopes using the expression f{v,N) = {- — 5—^)^ and 

(ii) /(w,iV)'s hence found are plotted linearly with 1/N. 

Therefore for each v, a functional form of f{v,N) is obtained in terms of 1/iV. Then for each length, the slopes 
and intercepts from the functional forms of f(v,N) for different v are plotted as a function of v to obtain the final 
fL{v,N) expression. For the log interaction the form is found to be f^°^{v,N) = [2?;oz(l + §7^)] where aoz ^ 16 for 
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(b). Ln Z- L Ln X Vs Ln L (Small N) 
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FIG. 1: LnZl{N) - LlnX is plotted with InL for the log interaction and oz models for (a). N = 1000000 and (b). = 1. 



L — 160 and for the oz model it is f°^{v, N) — [^'02(1 + ^)] • This results in two observations: (i). The base pairing 
interaction strength in 0] is doubled. This implies, v in the log interaction model is re-scaled to twice the v in 
the oz model which is similar to the way in which the parameter a in the linear interaction model re-scaled v to 
v/{l — a)^. (ii). An additional dependence of the combinatoric factor on N is found (va/N) for all small and large 
values of N for the log interaction and oz models (which comes out to be same). Such a dependence has not been 
found before in these matrix models. 



Power law exponents 



In order to extract the exponents corresponding to secondary and tertiary contributions, the InL dependence in 
the partition function is studied. This analysis is done for L = 160 and is divided into small N {N = 1) and large N 
(iV = 100000) regions. 



Large N 



To study the secondary exponent, {InZ — LlnX) is plotted with InL (Figl(a)). (i) The plot is a straight line with 
slope —1.37 which approaches —1.5 when larger and larger values of N are considered. This implies at large N, the 
major contribution is from the genus zero (secondary) structures and the power law L~^^^ holds true as for the oz 
and linear interaction models [3|-|6[. (ii) For v = 1, the oz and log interaction curves are distinctly different. On 
substituting v = 1/2 in the log partition function, the two model curves coincide. So for large N and v — 1/2, the 
structures in the log interaction model reduce exactly to the structures found in the oz model. Thus the effect of 
interaction at large N is visible in terms of v only, (iii) The slopes of the curves for different w's are nearly —1.37 
implying that the slope (i.e. the exponent —3/2) is independent of v. (b). Next, to study the tertiary exponent, 
{InZ — LlnX + {3/2)lnL) is plotted with L which shows: (i) the structure of the curves for the two models is exactly 
the same. So the genus contribution in the large N limit has the same form for both the models, (ii) The extrapolation 
of curves on y axis gives the value of fcg (for large N, see [3,]). For v = 1, these are different for the two models showing 
that the coefficients kg's will be different. For large values of v, the points for the even and odd lengths (at small 
lengths) split up for both the models (also seen at small N). The splitting is more for log model and may be due to 
different fcg's which seem to depend upon v. 



Small N 



{LnZ — LlnX) is plotted with InL (Fig 1(b)) and the following observations are made: (i) The effect of interaction 
in the small N region is due to v and N where N contributes dominantly while the contribution of v is very little, (ii) 
The different v curves for the two models clearly indicate that no value of v will ever reduce the log model to oz (or 
vice versa). This is a vital difference between the two models to establish their uniqueness particularly at small N. 
(iii) The plot is no longer a straight line with slope —1.5 but a U shaped curve which also includes the contribution 
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from crossed structures, (iv) The length at which secondary contribution becomes less dominant than tertiary (given 
by minima of the curve) is larger for the log interaction model. Therefore the effect of interaction is mainly on the 
secondary structures for a given length ll| . The large number of structures obtained from the asymptotic analysis 
of the partition function can be attributed to the increased structure space of all the n-punctured Riemann surfaces 



SPECIFIC HEAT 



The specific heat is defined as C(^y^i^^ = — T(^^) where F is the total free energy of the polymer chain for a 
given length L. The calculations (in Fig. 2) are performed for the log interaction and oz models for L — 160 and 
iV = (1, 3, 6, 10, 100000) with v = e~'^l^'^ (where Boltzmann constant and base specific binding energy ratio e/fc = 1). 
The following observations are made: (i) with T: For a fixed length, the peak value is maximum for = 1 and 
decreases (almost half) as larger values of N are considered for the two models. The peak decreases upto certain 
value of N after which it becomes nearly constant however large is the value of iV for both the models, (ii) dCy/dT 
with T: The first derivative of specific heat shows a kink for small N (for both the models) whereas for large iV's, no 
such behavior is seen. Further, the kink shifts to smaller temperatures as length is increased, (iii) (PCv/dT^ with T: 
At small T's, the curve shows (unusual) double peaks for = 1 for the log interaction model. There is a systematic 
conversion of the double peak into becoming a single peak as N is increased slowly. For large N, there is a peculiar 
kink present in the lower part of the curve for both the models. The kink is slightly more pronounced in the log 
interaction model than the oz. A similar such kink is visible (for largest length considered L = 1024) in the (PCy/dT^ 
verses T curves of the model in 7] (Pagnani et al) which discusses a disordered (glassy) statistical model of RNA 
secondary structures. The specific heat analysis visibly presents the peculiar differences in the characteristics of the 
two models at small and large A^'s with = 1 showing an unusual double peaked behavior. 



TAPPING 



In tapping, the matrix is coupled to an external source, which is then removed and the number of different config- 
urations (multiple solutions) are counted (Deo Q). The limit of external source — > gives different moments which 
may result in different partition functions and hence free energies (considering different tappings explores entire space 
of the configurations). For the log model, the form of potential is V{a) ~ [f — ^(logcr^)] with as the eigenvalues 
of fj. Introducing a linear matrix source a /y in the action of Eq. (1) (after solving), the saddle points are found from 

the equation " — = to be 



Ml/ V y A* 



1 I 4t"_ 

(15) 



2 

where ix=\/v and t" = 1/w. Depending upon small or large 1/y, the saddle points are 




- + i( — );->0. (16) 
^ 2 y 



^ = + ); -<o. (17) 

" ^ 2 y 



The action S{a) = §ct^ — ^log{a'^) — -a therefore becomes 



f" f" h" 1 

S[at) = -[l-log{-)]T\ — . (18) 
2 ^ y ^ y 



Thus the set of moments grows exponentially as 2^/^. This suggests the intriguing possibility that log interaction 
matrix model (and also oz) represent a class of 'glassy' matrix models. 
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FIG. 2: The figure shows specific heat Cv, C{, (inset) and C" as a function of T for the logarithmic interaction (solid line) and 
OZ models (dashed line) for different values of A^. 



CONCLUSIONS 



The Letter studies a Penner type logarithmic interaction in the framework of a random matrix model of RNA 
folding and its structure combinatorics. This is the first instance where Penner matrix models have been applied to 
the study of RNAs. The asymptotic analysis suggests: (i) There exists different regions of N in which the structural 
properties of the log interaction model (and also the oz) are different, (ii) A dependence of the combinatoric factor 
on the matrix size N is found for the matrix models of RNA with and without interactions which has not been found 
before, (iii) The effect of interaction is visible in v and N which is given by the re-scaling of v to twice that in 
the OZ (for large N region) and dominated by N (in the small N region). The substitution v/2 in the log partition 
function at large N, reduces the structures to that found for the oz model, (iv) The Lr'S-^/'^ behavior for large N 
in the log interaction model is in good agreement with the oz model. For small iV, these exponents change as the 
secondary region is stretched longer while the tertiary region gets shortened. The effect of interaction is thus visible 
in the way the structures are distributed for a given length which may represent a new universality class, Fig. 1(b) 
0. The specific heat clearly highlights the different structural behavior in the two regions of N with an unusual 
double peak for iV = 1 and a kink for large N . The tapping explicitly shows the presence of multiple solutions in 
these matrix models. An important necessary direction in these models is the derivation of the generating function 
with any N't' . This will help clarify some unresolved issues such as the genus characterization (which at the moment 
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is indistinguishable from effect of interaction both of which are contained in A''). The matrix model of RNA with 
logarithmic interaction provides a larger class of interacting RNA structures than the matrix model of [^-Q . This is 
because an additional (and therefore complete) dependence of the structures on n (other than vertices and edges) in 
the model will be possible in terms of all the n punctures of the Riemann surfaces. 
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